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Abstract 



^ \ The phase structure of lattice QCD with two flavors and Wilson fermions is studied an- 

O ' alytically. At /? = we obtain rigorous lower and upper bounds for the critical hopping 

^-i • parameter kc{^) from a convergent hopping parameter expansion to infinite order. The 

t~^ ! result supports the value fcc(O) = | observed in Monte Carlo simulations. 
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1 Introduction 

The outstanding properties of QCD are asymptotic freedom and confinement. Although 
a proof that QCD is confining does not exist, there are strong indications that QCD has 
this property and we assume it to be true in the following. The coexistence of asymptotic 
freedom and confinement means that the effective quark-gluon and gluon-gluon coupling 
changes continuously from weak coupling at short distance to strong coupling at long dis- 
tance. At long distance quarks and gluons are confined inside color singlet hadrons. All 
present known strong interaction phenomena support the confining picture. In addition 
hadrons built from the two lightest quarks exhibit an approximate flavour symmetry which 
is spontaneously broken. Due to the confinement, a perturbative description of these phe- 
nomena is not possible. The description of the long distance strong forces requires non- 
perturbative methods. Lattice QCD formulates the theory non-perturbatively and allows 
us to test confinement and to determine numerically properties of low energy hadrons. 

Some useful qualitative informations in lattice QCD can be obtained in the infinitely 
strong bare gauge coupling limit /3 = 0. At infinite gauge coupling confinement and spon- 
taneous breakdown of chiral flavour symmetry are manifest. The spontaneous chiral sym- 
metry breakdown and the emergence of massless Goldstone bosons explain the lightness of 
the pion. The small but non vanishing mass of the pion is due to the non exactness of the 
flavour symmetry due to the small but finite light quark masses. 

In the past many Monte Carlo simulations of lattice QCD with Wilson fermions has 
shown the existence of a critical line fcc(/3) in the k — (5 plane, where the mass of the Gold- 
stone boson vanishes. For any given coupling (3 there exists a critical hopping parameter 
kc{f3). The continuum limit is approached along this line which ends at kc{oo) = i. In the 
strong coupling regime, at /3 = 0, Monte Carlo simulations indicate that kc{0) = 4. Some 
strong coupling expansion analysis give also some good argument to support this feature 

Up to now there is not a rigorous proof that kc{0) = j for Wilson fermions. On the 
other hand, at infinite gauge coupling an important simplifiction of the system takes place 
which makes feasible a convergent hopping parameter expansion analysis to infinite order. 
In fact, at /? = all gauge variables on different links are independent and can be integrated 
out explicitly, leaving us with a pure fermion problem which is however still not trivial. 

In this work we analyse rigorously the behaviour of the Goldstone boson spectrum at /3 = 
with the method of the hopping parameter expansion. Our analysis is based on a random 
walk representation of the meson propagators |^]. In particular this allows us a rigorous 
estimation of the mass of the pseudoscalar isotriplet and of kc{0). A full hopping parameter 
expansion to infinite order of perturbation is a very difficult problem because there are no 
known general methods to identify all diagrammatic contributions to the expansion for a 
given order of the perturbation expansion. However, there is an easy way to determine 
an upper and a lower bound to kc{0) to infinite order of the perturbation expansion by 
neglecting and overcounting, respectively, a class of diagrammatic contributions of the 
expansion which are very difficult to be identified. This simplifiction leads to an upper 



and a lower estimation of the mass of the Goldstone boson (the pseudoscalar isotriplet) 
expressed simply by geometric series of the hopping parameter k. These bounds allow us 
to determine bounds on the critical hopping parameter kc{0) from the condition that this 
mass has to vanish. 

2 The Wilson lattice action 

There is a SU(3) color matrix U(b) in the fundamental representation defined on each 
oriented lattice bond b. Our convention is that 

U{-b) = U\b). (1) 

An oriented path uj on the lattice is a set of bonds 

uj = biUb2U ...Ubn (2) 

such that the endpoint of 6j is the start point of 6j+i for 1 < i < n — 1. We can associate a 
SU(3) color matrix with m by defining the path ordered product 

U{i0) = U{bi)U{b2)...U{bn) (3) 

The spin matrices are defined in terms of 7-matrices by 

r(h\ — J T^ = r + jf^ if b in +fi direction , . 

I f '^ = r — 7^^ if b in —fi direction 

The 7-matrices are hermitian 4x4 matrices, satisfying {7'^,7^} = 25^'^. 
The Wilson action is then defined on a lattice A by 

S = l3j2TrU{dp)-k J2 V'(2;)r(6)f/(6)V'(y) - ^ V'WV'W (5) 

pCA b=(xy) x<=A 

where p represents a plaquette. The quark fields are represented by the anti-commuting 
variables ^{x) and 'il>{x) which transform under the 3 and 3 representation of color. We 
consider only two flavors. 



3 The meson propagator 

In order to evaluate masses of mesons we have to consider the propagator of these states. 
The masses are obtained by studying the long distance behaviour of these propagators. A 
propagator of a meson is given by 

{M{x)^M{y)) = ^ f[dtP][di;][dU]M{x)^M{y)e-^'^^''^''^^ (6) 

where Z is the partition functional and M(x) = ip{x)M'ip{x) represents a meson operator 
(M is some matrix in spin, flavour and color space). The hopping parameter expansion 



is obtained by breaking the action S into two parts: the unperturbated part 5*0 and the 
perturbation Sj. 

So = -Y^ ip{x)'il^{x) 

Si = -k Y. i;{x)T{b)U{b)i^{x). (7) 

b={xy) 

Here the plaquette term is not present because we are at /? = 0. Expanding the exponential 
in eq. (6) in term of the perturbation and integrating out the gauge degrees of freedom U 
we obtain an expression for the the meson propagator 



(M(x)tM(y)) 



a;"'':xi— >y 



= Y CMxA{iv+Uiv-) (8) 

where Cm is some overah constant amphtude which does not affect the meson mass, u;"*" 
is a path on the lattice, uj~ = —uj~^ is a path in the opposite direction of u;"*" and L is the 
number of closed loops formed by the path u;+ U uj~ . Here T{uj) and U{uj) represent the 
ordered products of r(6) and U{b) on the bonds b G u, respectively. We do not give a proof 
of (8) since it is standard @, |J. 

4 Meson masses from the hopping parameter expansion 

In order to compute the meson masses we consider the static propagator 

G*^(t) = X:(M((0,f))tM((t,x))). (9) 

X 

The lowest order diagram representing a static propagator is a double fermion path uiq in 
the time direction (Fig.l). Since there is a non-zero probability for a transition from the 
lowest order state to a more complicate state, the full propagator is given by a sequence 
of exited states connected by lowest order states (Fig. 2). Each double path uj ^ ujq can be 
viewed as a space-time process contributing to the excitation. These excitations renormalize 
the mass of the static propagator. 

The unrenormalized mass of a meson can be found from the properties of the lowest order 
diagram of the static propagator which is given by 

G^'it) = Cm X exp(-mot) = Cm x {4ky. (10) 

for r = 1 in (4). This tells us that the unrenormalized mass is mo = —log{'ik'^). Having 
computed the unrenormalized static propagator, we will now compute the renormalization 
produced by intermediate exited states. Fig. 3 shows us such an excitation, with initial 
point w and final point z, where the static unrenormalized propagator arrives and departs. 



For fixed w and z, we will sum over all intermediate exited states to obtain a total weight 
for the event. We denote this weight by D [w, z). The full propagator takes the form Q 

(M(x)tM(y)) = Yl Gf{x,wi)D^'{wi,zi)Gf{zi,W2)D^'{w2,Z2)x 

Wi,Zi(i=l...n) 

X...xD^'{Wn,Zn)G^\zn,y) (11) 

where the points Wi and Zi are required to lie between x and y and to be ordered so that 
each Zi is "later" than Wi, which is itself "later" than Zj-i. This last equation represents 
the picture in Fig. 2. 

In order to find the renormalized mass of the propagator, we consider the full static prop- 
agator 

G''{t) = ^(M(x)tM(y)) 



x,y 



(12) 
yo-xo=t 

where xq and yo represent the time component of the lattice points x and y, respectively. 
Similarly we defineg 



x,y 



(13) 
yo-xo=t 

Here we have normalized the excitation relative to a static unrenormalized propagator over 
the time interval t. From this definition and eq. (11) we obtain 



G'""' (t) 



g-moi J2 D^Hh - 


-Si). 


.D'^it 


Si,U{i=l...n) 






g-moi gPAft 







(14) 



where ij and Sj are the corresponding time coordinates of Zi and Wi. The exponential 
representation of (14) is up to irrelevant boundary effect exact. So the renormalized mass 

niM is 

ruM = mQ-pM- (15) 

The leading order oi pM can just be written in the form 

m= E {G^\t'))-^D'\Q,z)= Y: {^k')-''D'^iO,z) (16) 

z={t'>0,z) z={t'>0,z) 

To proof this last equation we have to insert it into the last term of (14) and to expand the 
exponential function containing pM , then we have to compare the result with the first term 
of (14). The excitation term D [w, z) is defined starting from eq. (8) in the following way: 

D''{w,z)= E 4^ (17) 

where the sum is over all closed paths to = lo^ U uj^ from w to z and return. A{uj) denotes 
the amplitude of an intermediate excitation diagram and Cm is the overall constant of eq. 
(8) and (10). 

^Notice that {G^' {t))-^ = c^e'"«* = -^{Ak")-' 



5 Goldstone boson mass 

The critical hopping parameter kc{0) can be evaluated from the condition that the mass 
of the Goldstone boson (the pseudoscalar isotriplet) vanishes. The mass of the Goldstone 
boson is expressed in term of a hopping parameter expansion serie by eq. (15-17). The 
Goldstone boson is characterised by the pseudoscalar isotriplet bound state matrix M = 
Mcb = 75^ for I = 1,2, 3, where T are the isospin SU(2) generators. 
To evaluate the correction to the unrenormalized mass mo due to the excitation pM we 
have to identify all closed paths cv from the point to the point z and to sum over all 
points z = {t > 0,z) with positive time in eq. (16). An excitation is a closed path. At 
/? = we do not have plaquettes in the expansion because the plaquette expectation value 
vanishes, therefore an excitation is a closed path consisting of a double line uj = to^ U uj~ 
with L0~ = —uj^ composed by a set of connected bonds 



LO 



6J+ UUJ' 



u ^^ 

bieu}+ 



U 



U (-b^) 

bieu)+ 



(18) 



An excitation has to be irreducible, which means that cutting a double bond 6 U — 6 of the 
path LO in the time direction, the resulting two paths have some time overlap. Reducible 
excitations can be split into several time separated excitations, therefore to avoid over- 
counting we have to consider only irreducible excitations. 

Using the following two properties for r = 1 of the gamma matrices r(6)r(6) = 2r(6) and 
r(6)r(— 6) = the reader can easily convince himself that the amplitude A{uj) of an exci- 
tation of the Goldstone boson can be expressed by a product of terms depending on the 
order of the bonds which compose the path to 

A{uo = to+ U L0-) = Cmg, X Ao X [] F{bi, b,+i)k^ (19) 

biCui+ 

where Aq and Cmoh ^^^ some constants and 

if bi+i = -hi 
F{h,hi+i) = { 2 iihi+i^h (20) 

if fej+i II hi 

The hopping parameters is squared in (19) because the product is only on the half path ui^ 
and all the bonds of lo^ = —lo^ are paired with bonds in lo^ . 

There are two classes of irreducible paths lo. The first contains paths with all bonds bi 
perpendicular to the time direction. These paths are trivially irreducible, because they 
have no bonds in the time direction. It is very easy to identify and count them. The second 
class contains irreducible paths which have some bond parallel to the time direction (Fig. 
4). These paths are very difficult to be counted because they are irreducible and therefore 
they can not be decomposed into a product of spatial excitations separated in time and 
time-like double bonds. 

Because the function F in (20) is positive we can however undercount and overcount the 
irreducible paths contributing to the sum in (17), by neglecting all paths of the second class 
and, respectively, considering all possible paths (reducible and irreducible) with bonds in 



all directiond^. We obtain in this way an underestimation and an overestimation of PMch ' 

pXT = ^oxJ] W C{hA+i)F{hA+i)k^ (21) 

S'cShi£S' 

PmZ = Ao>^Y.T. n C(6„6,+i)F(6„6,+i)fc2 (22) 

S'CST'cTbiGS'UT' 

where the set S represents all bonds perpendicular to the time direction and the set T all 
bonds parallel to the time direction. The term (4/c^)~* in the sum (16) is not present in 
(21) because there is no time extension in the sum in (21). In (22) it is implicitly present in 
the overestimated form (A;^)~* . The mapping C(6j, &i+i) controls if the sequence of bonds 
forms a connected path 

„,, , \ _ J if ^«+i ^^'^ ^i ^^6 sequentially disconnected , . 

[1 if 6j+i and 6j are sequentially connected 

We define two bonds bi and ftj+i to be sequentially connected if the end point of 6j is the 
start point of 6i+i. It is easy to see that the set of all possible paths which contribute 
non-trivially to eq. (21-22) forms a 3-dimensional and, respectively, a 4-dimensional tree. 
The number of branches of these trees counts the number of different paths. For a n- 
dimensional tree each bond can be chosen perpendicular or parallel or antiparallel to the 
previous one. For all these possibilities the function F is known and given by (20). Because 
of the particular form of F we obtain for (21) and (22) two geometrical serieqj: 

oo 

(24) 
(25) 

n=l 

where d=4 is the dimension of the space-time lattice and Aq = ^ is evaluated explicitly 
from eq. (8). 

The final result is an over and upper bound on the Goldstone boson mass. From eq. (15) 
we obtain 



under _ 
PMob - 


= Aoyx^d-i)er 

n=l 


^over 


oo 

= Ao)^iMk'r 



^ ^ < ™„„.(*. ,3 = 0) < - lo.(4.^) - ^^^),, 



los("') - 7— TTU S ™«o.(*. /3 = 0) < - 106(4*") - , \,, ' .„,, (26) 



From the condition inMob = we obtain two bounds on kc{(i = 0) 

0.231 < kciP = 0) < 0.264. (27) 

This rigorous result supports the observed value in lattice Monte Carlo simulations. 

The same method can be used for finding an upper and a lower bound on the masses of the 



^We recall that considering in the sum (17) a reducible paths has the effect of overcounting some irre- 
ducible paths. 

^We consider, for example, eq. (22). For a given bond bi we can choose &i+i perpendicular (2(d — 1) 
possibilities) or parallel (1 possibility) or antiparallel (1 possibility) to bi. Using the value of the function F 
we obtain for these two bonds: 2 x 2(d — 1) +4x1+0x1 = 4d. Eq. (21) is analogous to eq. (22). 



p and u) mesons. The p and u mesons are characterised by the matrices Mp = ^^T^ (j = 
1, 2, 3) and M^ = 7-^ (j = 1, 2, 3), respectively. One has to find the corresponding functions 
F in eq. (20). For other mesons the method can not be appUed because the function F 
vanishes. 
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Figure Caption 

1. A static unrenormalized propagator. 

2. An excitation of the static unrenormalized propagator. 

3. A close up of an intermediate excitation. 

4. A typical irreducible excitation with some bond in the time direction. 
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